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■ Abstract 

A global fit of parameters allows us to pin down the Hidden Local Symmetry (HLS) 
effective Lagrangian, which we apply for the prediction of the leading hadronic vacuum 
polarization contribution to the muon g — 2. The latter is dominated by the annihilation 
channel e^e~ — )• tt'^tt", for which data are available by scan (CMD-2 & SND) and ISR 
(KLOE-2008, KLOE-2010 & BaBar) experiments. It is well known that the different 
data sets are not in satisfactory agreement. In fact it is possible to fix the model parameters 
without using the vr+vr^ data, by using instead the dipion spectra measured in the r-decays 
together with experimental spectra for the t:^^, tjj, tt^tt^tt^, K^K~, K^K^ final states 
supplemented by specific meson decay properties. Among these, the accepted decay width 
for — )• e+e^ and the partial widths and phase information for the uj/cl) — )• vr+vr" 
transitions, are considered. It is then shown that, relying on this global data set, the HLS 
model, appropriately broken, allows to predict accurately the pion form factor up to 1 .05 
GeV It is shown that the data samples provided by CMD-2, SND and KLOE-2010 behave 
consistently with each other and with the other considered data. Consistency problems 
with the KLOE-2008 and BaBar data samples are substantiated. "All data" global fits are 
investigated by applying a reweighting to the conflicting data sets. Constraining to our 



best fit, the broken HLS model yields a*f = (11 659 169.55 + 



+1.26 
-0.59 



± 5.2lth) 10 



-10 



associated with a very good global fit probability. Correspondingly, we find that Aa^ 
o^^P — a*J^ exhibits a significance ranging between 4.7 and 4.9 a. 
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1 Introduction 



The theoretical value for muon anomalous magnetic moment is an important window 
in the quest for new phenomena in particle physics. The predicted value is the sum of several 
contributions and the most prominent ones are already derived from the Standard Model with 
very high accuracies. The QED contribution is thus estimated with an accuracy of a few 10^^ 

121 [3]| and the precision of the electroweak contribution is now of order 10^^^ H. The 
light-by-light contribution to is currently known with an accepted accuracy of 2.6 x 10~^° 
ISJ. 

Presently, the uncertainty of the Standard Model prediction for is driven by the uncer- 
tainty on the leading order (LO) hadronic vacuum polarization (HVP) up to ~ 2 GeV (HEII- 
This region is covered by the non-perturbative regime of QCD and the leading order HVP 
{LO — HVP) is evaluated by means of : 



which relates the hadronic intermediate state contributions {Hi, i = 1 ■■ - n} to the annihilation 
cross sections cr(e+e^ — )• Hj) = auX^)- ^i^) ^ known kernel [4J enhancing the weight of 
the threshold region s^^ and Scut is some energy squared where perturbative QCD starts to 
be applicable. In the region where perturbative QCD holdfl its contribution to carries an 
uncertainty of the order of a few 10~^^. 

Up to very recently, the single method used to get the a^{Hi) 's was to plug the experimental 
cross sections into Eq. ([T]). Among the most recent studies based on this method, let us quote 
HI IS 121 HJ. When several data sets cover the same cross section crj/. (s), Eq. is used 
with some appropriate weighting of the various spectra, allowing to improve the corresponding 

On the other hand, it is now widely accepted that the Vector Meson Dominance (VMD) 
concept applies to low energy physics [fTOlfTTIl . VMD based Effective Lagrangians have been 
proposed like the Resonance Chiral Perturbation Theory or the Hidden Local Symmetry (HLS) 
Model; it was soon proved [12J that these were equivalent. Intrinsically, this means that there 
exist physics correlations between the various e+e^ — )■ Hj annihilation channels. Therefore, it 
becomes conceptually founded to expect improving each a^(ifj) by means of the data covering 
the other channels e+e^ — )■ Hj (j ^ %). 

This is basically the idea proposed in [T3l relying on the HLS model fT4l [151. Using a 
symmetry breaking mechanism based on the simple BKY idea [|T6 l and a vector meson mixing 
scheme, the model has been developed stepwise [fTTl [T8l [T9l |20l and its most recent form |fT3l 
has been shown to provide a successful simultaneous description of the e+e^ annihilation into 
the TT+TT^, 7r°7, 777, tt+tt^tt^, K^K~ , K^K^ final states as well as the — )■ z/^tt^tt^ decay 
spectrum. Some more decays of the fonr@ V — t- P7 or P — > 77 are considered. 



'The charmonium and bottomium regions carry uncertainties also in the range of a few 10 

^We denote by y or P resp. any meson belonging to the (basic) vector or pseudoscalar lowest mass nonets. 
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As higher mass meson nonets are absent from the standard HLS model, its energy scope is 
a priori limited upwards by the (f) meson mass region (~ 1.05 GeV). However, as this region 
contributes more than 80 % to the total HVP, improvements which can follow from the broken 
HLS model are certainly valuableEl. 

The global simultaneous fit of the data corresponding to the channels quoted above allows 
to reconstruct the various cross sections (T//. (s) taking automatically into account the physics 
correlations inside the set T-L = {Hi} of possible final states and decay processes. The fit param- 
eter values and the parameter error covariance matrix summarize optimally the full knowledge 
of "H. This has two important consequences : 

• One should get the {a^(Hi),i = 1, - ■ - n} with improved uncertainties by integrating the 
model cross sections instead of the measured ones. Indeed, the functional correlations 
among the various cross sections turn out to provide (much) larger statistics in each 
channel and thus yield improved uncertainties for each a^{H.i). 

• When several data samples cover the same process Hi, one has a handle to motivatedly 
examine the behavior of each within the global fit context. Stated otherwise, the issue of 
the consistency of each data set with all the others can be adressed with the (global) fit 
probability as a tool to detect data samples carrying problematic properties. 

Up to now, the broken HLS model (BHLS) [fT3l - basically an empty shell - has been 
fed with all existing data set^ for what concerns the annihilation channels 7r°7, 777, tt+tt^tt^, 
K+K-, K^Tf, with the spectra from ALEPH CLEO Ell and BELLE ^ for the r 
dipion deca)0 and with the VP'y / P77 partial width information extracted from the Review 
of Particle Properties (RPP) [|24|. This already represents more than 40 data sets collected by 
different groups with different detectors; one may thus consider that the systematics affecting 
these data sets wash out to a large extent within a global fit framework. 

For what concerns the crucial process e+e^ — t- tt+tt^, the analysis in [fT3l only deals with 
the data sets collected in the scan experiments performed at Novosibirsk and referred to glob- 
ally hereafter as NSK Il25ll26ll27ll28ll29l . The main reason was, at this step, to avoid discussing 
the reported tension [[30l [H between the various existing tt+tt" data sets : the scan data sets 
just quoted, and the data sets collected using the Initial State Radiation (ISR) method by KLOE 
ll3Tl[32l and BaBar [|33l[34ll . not to mention the pion form factor data collected in the spacelike 
region (351 [31. 

It has thus been shown that the global fit excluding the ISR data sets, allows to yield a 
splendid fit quality; this proves that the whole collection of data sets considered in lfT3l is self- 
consistent and may provide a safe reference, i.e. a benchmark, to examine the behavior of other 
data samples. 

Using the fit results, the uncertainty on the contribution to a^ of each of the annihilation 
channels considered was improved by - at least - a factor of 2, compared to the standard 

^The broken HLS model does not include the An, 5tt, 6tt, t^tttt and ujtt annihilation channels. Therefore, the 
(small) contribution of these missing channels [13 ] to should be still evaluated by direct integration of the 
experimental cross sections; up to the mass, this amounts to [l 31 (1.55 ± 0.57tot) 10^^°. 

"^The full list of data sets can be found in |fT9l or llT3l together with a critical analysis of their individual 
behavior 

^The energy region used in the fits has been limited to the [2m7r, 1 GeV] interval where the three data sets are 
in accord with each other This should lessen the effect of some systematic effects. 
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estimation method based on the numerical integration of the measured cross sections. For 
the case of the tt+tt^ channel, the final uncertainty was even found slightly better than those 
obtained with the standard method by merging scan and ISR data, i.e. a statistics about 4 times 
larger in the tt+tt^ annihilation channel. 

The main purpose of the present study is an update of the work in [fT3l aiming at confronting 
all scan (NSK) and ISR (from BaBar and KLOE) - and even spacelike (3511361 - tt+tv- data 
and reexamine the reported issues ll30l [8l. The broken HLS model described in |[13l happens 
to provide a tool allowing to compare the behavior of any of these tt+tt^ data sets when con- 
fronted with the 7r°7, 777, tt+tt^tt^, K^K~ , K^K^ annihilation data and with the r dipion 
spectra. Indeed, these data alone, supplemented with some limited information extracted from 
the Review of Particle Propertied (RPP) fi24il . allow to predict the pion form factor with a 
surprisingly good precision. The additional RPP information is supposed to carry the Isospin 
Breaking (IB) information requested in order to derive reliably the tt+tt information from the 
knowledge of the ti^tx^ spectrum. 

We also take profit of the present work to update the numerical values for some contribu- 
tions to the muon anomalous moment a^, all gathered in Table 10 of [13J. Thus, we update 
the QED entry by using the recent spectacular progress by Aoyama, Hayakawa, Kinoshita and 
Nio (H 0. They have been able to perform a complete numerical calculation of the 5-loop 
QED corrections to Og and a^. On the other hand, the electroweak contribution, which depends 
on the Higgs mass at 2-loops is now better known if we accept that ATLAS [|37l and CMS [[381 
have observed at the LHC the Higgs boson at a mass of about 125 GeV in a narrow window. 
Using this information slightly changes the central value as well as the uncertainty of the EW 
entry. We also have reevaluated the higher order HVP contribution (HO) within the standard 
approach based on all tt+tt^ channels {i.e. all scan and ISR data). 

The paper is organized as follows. Section [2] reminds (from [T3l) a few basic topics con- 
cerning the TTTT channel description within the BHLS model. In Section [3l the detailed frame- 
work - named "r-i-PDG" - used to study the differential behavior of the scan and ISR data 
is presented. Thanks to the (wider than usual) energy range covered by the BaBar spectrum 
[|33ll34l , a detailed study of the tt+tt" spectrum in the region can be performed for the first 
time. This leads to update the -> tt+tt^ treatment within our computer code; this is empha- 
sized in Subsection l3.3[ In SectionlH one confronts the "r-i-PDG" predictions with the available 
scan (NSK) and ISR data samples; it is shown that the NSK data and both KLOE data samples 
(referred to hereafter as KLOE08 [31] and KLOEIO [32]) have similar properties while BaBar 
behaves differently, especially in the p — uj interference region. Section |5] reports on the global 
fits performed using the various tt+tt^ data samples each in isolation or combined. Subsection 
I5.1.4l is devoted to studying the region of the pion form factor. In Section |6l we focus on the 
consequences for the muon anomalous moment of the various scan and ISR tt+tt^ spectra 
and compare results with the BNL [|39ll40l measurement. The tt+tt^ intermediate state con- 
tribution to from the invariant mass region [0.630, 0.958] GeV is especially considered as 
it serves to examine the outcome of various fits with respect to the experimental expectations. 
Finally, Section|7]is devoted to conclusions. 

^We occasionally refer to the RPP as Particle Data Group (PDG). 
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2 A Brief Reminder Of Concern for the tttt Channel 



2.1 The Broken HLS Lagrangian 

The general framework of our analysis is the resonance Lagrangian extension [[TOl [T2| of 
Chiral Perturbation Theory (ChPT), the latter being the systematic and unambiguous approach 
to low energy effective QCD [|4T1 142]| . This approach can also be viewed as an attempt to 
implement the VMD model in a way consistent with the chiral structure of QCD. Like in ChPT, 
the basic fields are the unitary matrix fields ^l^h = exp [±i P/ fjr], where P = Pg, + Pq is the 
SU (3) matrix of pseudoscalar fields, with Pq and Pg being respectively the basic singlet and 
octet field matrices. The Hidden Local Symmetry (HLS) ansatz [T4l[T5l is an extension of the 
ChPT non-linear sigma model to a non-linear chiral Lagrangian based on the symmetry pattern 
Ggiobai/i^iocai, where G = SU{3)l ® SU{3)r is the chiral group of QCD and H = SU{3)v 
is the vector subgroup. The hidden local SU{3)v requires the spin 1 vector meson fields, 
represented by the SU{3) matrix field V^, to be gauge fields. The corresponding covariant 
derivative reads = — igVj^i and can be naturally extended [[TSl in order to include the 
couplings to the electroweak gauge fields A^, and W^^. As the coupling to the Z boson is 
not relevant in our context, it is omitted in the following. The unbroken HLS Lagrangian is 
then given by £hls = + ajCy, where 



with L = [Dfj_C,L]^l ^iid ^ = [Dfj.^R]^R'^ a is a basic HLS parameter not fixed by the theory 
which should be derived from data. From standard VMD models, one expects a ^ 2. 

In fact the global chiral symmetry Ggiobai is well known not to be realized as an exact 
symmetry in nature which implies that the ideal HLS symmetry is evidently not a symmetry 
of nature either. Therefore, it has obviously to be broken appropriately in order to provide a 
realistic low energy effective theory mimicking low energy effective QCD. 

Unlike in ChPT where one is performing a systematic low energy expansion in low mo- 
menta and the quark masses, here one introduces symmetry breaking as phenomenological 
parameters to be fixed from appropriate data. Indeed, a systematic low energy expansion a 
la ChPT ceases to converge at energies above ~ 400 MeV, while we attempt to model phe- 
nomenology up to, and including, the (f) resonance region. 

In our approach, the Lagrangian pieces in Eqs ^ undergo a first breaking mechanism 
named BKY, originally defined in [|16l . It has been reformulated in [|17il in order to avoid 
some undesirable properties. In its original form, this BKY breaking scheme covers only the 
breaking of the SU(3) symmetry; following [43J, the BKY breaking has been extended in order 
to include isospin symmetry breaking effects. This turns out to modify Eqs. Q by introducing 
two constant diagonal matrices Xa/v '■ 



and the (non-zero) entries in X^/v are fixed from fit to the data. The final broken HLS La- 
grangian can be written : 



^A/V 




(2) 




(3) 



HLS — 




(4) 
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where C'tHooft contains determinant terms L44il breaking the nonet symmetry in the pseu- 
doscalar sector. C'jj^g can be found expanded in fl3\. 

However, in order to account successfully for the largest possible set of data, isospin sym- 
metry breaking should be completed by the kaon loop mixing of the neutral vector mesons (p% 
LOj and (pi). This implies a change of fields to be performed in the Cj^ig Lagrangian. This 
mixing mechanism is briefly reminded just below. 



2.2 Mixing of Neutral Vector Mesons Through Kaon Loops 

It has been shown [[T9ll20l[T3Tl that £^^5 is insufficient in order to get a good simultaneous 
account of the e+e~ — )• tt+tt" annihilation data and of the dipion spectrum in the r"^ — )■ v^n^n'^ 
decay. A consistent solution to this problem is provided by the vector field mixing mechanism 
first introduced in [|T8l . 

Basically, the vector field mixing is motivated by the one-loop corrections to the vector 
field squared mass matrix. These are generated by the following term of the broken HLS 
Lagrangiar|Z|£^^5 : 



lag 

4:ZA 



I [p? + V2zv 0/] K-dK++ [p? -UJI + V2zv 0/] d i^"} , (5) 



where g is the universal vector coupling and the subscript / indicates ideal vector fields. 

Therefore, the vector meson squared mass matrix Mg , which is diagonal at tree level, un- 
dergoes corrections at one-loop. The perturbation matrix 5M'^{s) fTH [191 120'| depends on the 
square of the momentum flowing through the vector meson lines. The diagonal entries acquire 
self-mass corrections - noticeably the p° entry absorbs the pion loop - but non-diagonal entries 
appear and correspond to transitions among the ideal cu and meson fields which originally 
enter the HLS Lagrangian : li^^^s), Yip^^^s) and Iip^{s). These are linear combinations of the 
neutral and charged kaon loopq^- 

Therefore, at one-loop order, the ideal vector fields which enter C^jj^g are no longer mass 
eigenstates; the physical vector fields are then (re)defined as the eigenvectors of = Mq + 
5M^(s). This change of fields should be propagated into the whole broken HLS Lagrangian 
C'uj^g, extended in order to include the anomalous couplings |f45l as done in [T3]. For brevity, 
this Lagrangian is referred to as BHLS. 



2.3 The Vtitt and V - -f/W^ Couplings 

As the present study focuses on e+e^ — t- tt+tt^ data, it is worth to briefly remind a few 
relevant pieces of the C^jj^g Lagrangian. In terms of physical vector fields, i.e. the eigenstates 

' For clarity, we have dropped out the isospin breaking corrections generated by the BKY mechanism; the 
exact formula can be found in the Appendix A of IHJ. The parameter zv corresponds to breaking of the SU(3) 
symmetry in the Lagrangian piece Cy, while za is associated with the SU(3) breaking of Ca- zy has no really 
intuitive value, while za can be expressed in terms of the kaon and pion decay constants as za = [/if / f^Y- 

^Other contributions, like K*K loops, take place fT8',n[9| which are essentially real in the energy region up to 
the meson mass. These can be considered as numerically absorbed by the subtraction polynomials of the kaon 
loops. 
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of = Mq + 5M^(s), the Vnn Lagrangian piece writes : 



^(l+Si/) [{p° + [(1 - hv)Av - a{s)] u + /3{s) 0} • vr^ H 7r+ + |p- ■ 7r+ 5 7r° - p+ ■ vf- a vr^ 

(6) 

Sy and (1 — /ly) Ay are isospin breaking parameters generated by the BKY mechanism [[T3l . 
whereas a(s) and /3{s) are (complex) "angles" which define the physical neutral meson vector 
fields in terms of their ideal partners. Their expressions can be found in [JSJ . 
Another Lagrangian piece relevant for the present update is : 



■A 



-fpw W+-P- + W-- 



(7) 



where g2 is the weak coupling constant related with the Fermi constant [fT5B and Vud is the 
element of the (u, d) entry in the CKM matrix. The fv-yis) functions and fpw are the transition 
amplitudes of the physical vector mesons to resp. a photon and the W boson. At leading order 
in the breaking parameters, they are given by [|T3l : 



f> 



3 



^ , Ay a(s) V2zv ni \ 

1 + Sy + hy^ + + -^^^^^ 



1 + Sy + 3(1 - /iy)Ay - 3a(s) + V2zy 7(s) 



(8) 



-y2;zy + 3/3(s) +7(s) 



f,w ^ /; = agfl [1 + Sy] . 

Eqs. dll) and ([8]) exhibit an important property which should be noted. The functions 
providing the coupling of the physical uj and <\) mesons to a pion pair also come into ffry{s). 
Therefore, any change in the conditions used in order to account for the decays uj/cj) ^ tx^tx^ 
correspondingly affects the whole description of the e+e~ vr+vr" cross section. 



2.4 The Pion Form Factor 



We only remind the pion form factor in r decay and in e^t~ annihilation and refer the 
interested reader to [13] for detailed information. The pion form factor in the decay to 
Tx^Tx'^v^ can be written : 

1 



1--(1 + Sy) 



(9) 



where a and g are the basic HLS parameters [15] already encountered; Sy is one of the isospin 
breaking parameters introduced by the (extended) BKY breaking scheme. The other quantities 
are : 



Dpis) = s-ml-U'{s) 



(10) 



[ f; = agfl{l + J:v) 



ag'fl{l + J:v) 
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where Hwis) and Il'pp{s) are, respectively, the loop correction to the p"^ — transition 
amplitude and the charged p self-mass (see L13J). 

The pion form factor in e~^e~ annihilation is more complicated and writes : 



_ pe ( \ 9 pirn 



, (11) 



where (^pTrvr, gw-Kw and g^j.^, can be read off Eq. Q and the Fy are given by : 



Fv^is) = fv^{s) - Uv^is) , {V = p% ur, (pn) 



(12) 



with the fvyis) given by Eqs. ([8]) above and the ny^(s) being loop corrections [fT3l . Dp{s) 



npp(s) is the inverse p° propagator while D^{s) and D^{s) are the modified fixed 



width Breit-Wigner functions defined in [ 13]. 



3 Confronting the Various e+e — > tt+tt Data Sets 
3.1 The Issue 

Although the BHLS Lagrangian should be able to describe more complicated hadron pro- 
duction processes, in a first step one obviously has to focus on low multiplicity states, pri- 
marily two particle production but also the simplest three particle production channel e+e^ — > 
7r+7r^7r°. Four pion production, annihilation to KKn . . .are beyond the scope of the basic 
setup of the BHLS model. We expect that available data on the lowest multiplicity channels 
provide a consistent database which allows us to pin down all relevant parameters, such that 
our effective resonance Lagrangian is able to simultaneously fit all possible low multiplicity 
channels. In fact, what is considered are essentially all relevant annihilation channels up to the 
0; in this energy range, as already stated, the missing channels (47r, Svr, Gtt, r^vrvr, un) contribute 
less than 0.3% to a^'^. 

Our previous study [T3l has actually shown that the following groups of complementary 
data samples and/or RPP [24J accepted particle properties (mainly complementary branching 
fractions) support our global fit strategy : 

• (i) All e"'"e~ annihilation data into the 7r°7, 777, Tr+vr"??^, K+K-, K^K final states admit 
a consistent simultaneous fij^, 

• (ii) The r± VrTi^-n^ dipion spectra produced by ALEPH CLEO and BELLE 
(231, however limited to the energy region where they are in reasonable accord with each 
other (v/i < 1 GeV), 

• (iii) Some additional partial width from the P77 and VP'^ decays, which are independent 
of the annihilation channels listed just above, 

^ With - possibly - some minor tension between the Tr+Tr^Tr" data around the </) resonance and the dikaon data 
(see the discussions in |fT9l[T3]| ). 
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• (iv) Some information concerning the n^n^ decay, especially its accepted partial 
width r(0 — )■ Ti~Ti~^) ll24l . This piece of information is supposed to partly counterbal- 
ance the lack of spectrum for the e^e^ — t- tt+tt^ annihilation in the mass regior0. 

• (v) All the e+e^ — ^ tt+tt^ data sets (NSK) collecte(£l] by the scan experiments mounted 
at Novosibirsk, especially CMD-2 28] and SND ||29ll . 

They represent a complete reference collection of data samples and lead to fits which do 
not exhibit any visible tension between the BHLS model parametrization and the data (see for 
instance Table 3 in [13,|). This is worth being noted, as we are dealing with a large number of 
different data sets collected by different groups using different detectors and different accel- 
erators. The (statistical & systematic) error covariance matrices used within our fit procedure 
are cautiously constructed following closely the group claims and recommendations. There- 
fore, the study in [13J leads to think that the model correlations exhibited by BHLS reflect 
reasonably well the physics correlations expected to exist between the various channels. 

However, beside the (NSK) e+e^ — )■ tt+tt^ data sets collected in scan mode, there exists 
now data sets collected using the Initial State Radiation (ISR) method by the KLOE and BaBar 
Collaborations. All recent studies (see [301 [H, for instance) report upon some "tension" be- 
tween them. As this issue has important consequences concerning the estimate of the muon 
anomalous magnetic moment, it is worth examining if the origin of this tension can be identi- 
fied and, possibly, substantiated. Besides scan and ISR data, it is also interesting to reexamine 
ifTSl the pion form factor data collected in the spacelike region [[35l[36l within the BHLS frame- 
work; indeed, if valid, these data provide strong constraints on the threshold behavior of the 
pion form factor and, therefore, an improved information on the muon g — 2. 

3.2 The Analysis Method 

The BHLS model has many parameters and a global fit has to be guided by fitting those 
parameters to those channels to which they are the most sensitive. Obviously resonance pa- 
rameters of a given resonance have to be derived from a fit of the corresponding invariant mass 
region. Similarly, the anomalous type interaction responsible for 7r° — )■ 77 or the nnn final 
state ... are sensitive to very specific channels only. 

We also have to distinguish the gross features of the HLS model and the chiral symmetry 
breaking imposed to it. With this in mind, in our approach, the r decay spectra play a key 
role since the charged channel is much simpler than the neutral one where 7, to and are 
entangled by substantial mixing of the amplitudes, which are not directly observable. In the 
low energy region, below the kaon pair thresholds and the region, what comes into play is the 

form factor obtained from the r spectra. Together with the isospin breaking due to p° — w 
mixing - characterized by the branching fractions Br{uj — )■ 7r~^7[^) and Br(uj — t- e+e^), which 
in a first step can be taken from the RPP - the p"^ form factor should yield a good prediction 
for the e+e^ — > tt+tt^ channel. Data from the latter can then be used to refine the global fit. 
This will be our strategy in the following. 

'''The BaBar data [33J allow, for the first time, to make a motivated statement concerning how this piece of 
information should be dealt with inside the minimization code. This is discussed below. 
' ' The data sets [i25l collected by former detectors at Novosibirsk are also considered. 
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The annihilation channels referred to as (i) in the above Subsection as well as the decay 
information listed in (iii) have little to do with the e+e^ — ^ tt+tt^ annihilation channels, except 
for the physics correlations implied by the BHLS model. On the other hand, as long as one 
limits oneself to the region (2m^, 1 GeV), there is no noticeable contradiction between the 
various dipion spectra extracted from the — )■ i/^vr^vr^ decay by the various groups [1211 l22l 
|23J . Therefore, it is motivated to examine the behavior of each of the collected e+e^ — )■ tt+tt" 
data sets in isolation while keeping as reference the data corresponding to the channels listed in 
(i) - (iii). Stated otherwise, the data for the channels listed in (i) - (iii), together with the BHLS 
model, represent a benchmark, able to examine critically any e'^e^ — tt+tt^ data sample. 

On the other hand, the Isospin Breaking (IB) effects affecting the dipion spectra from r 
decay are independent of those in the e+e~ — )■ vr+vr^ annihilation and can easily be factored 
out. Therefore, it remains only to account for isospin breaking effects specific of the e+e^ — )• 
TT+TT^ channel, in a clearly identified way. 

A priori, IB effects specific of the e'^e^ — t- tt+tt^ annihilation are threefold and cover : 

• (j) Information on the decay p° — )• e+e^, 

• (jj) Information on the decay cu — )■ tt+tt^, 

• (jjj) Information on the decay — > tt+tt^. 

The importance of decay information on p° — )■ e+e^ to determine IB effects has been 
emphasized in only a few previous works ^ |20l [131. Within the BHLS model, the ratio 
fpy{s)/ fpw exhibits non-negligible IB effects for this particular coupling (see Fig. 1 1 in lfT3l ). 
They amount to several percents in a (threshold) region quite important for evaluating g — 2. 

There is certainly no piece of information in the data covered by the channels listed in (i) - 
(iii) above concerning the decay information (jj) or (jjj). In contrast, the vertex p°e+e~ is cer- 
tainly involved in all the annihilation channels considered. Imposing the RPP [24J information 
r(p° — )• e+e^) = 7.04 ± 0.06 keV is, nevertheless, legitimate because the channels (i) - (iii) 
do not significantly constrain the decay width p° — )■ e+e^. 

For the following discussion we define the branching ratio products F^^ = Bi{u — > e+e^) x 
Br{cu — )■ TT+TT^) and = Br(0 — )■ e+e") x Br(</) — )■ tt+tt"), which are relevant for the 
prediction of the e'^e" — )■ tt+tt^ channel. 

Two alternative analysis strategies can be followed : 

• (k) Use the so-called accepted values for the p° — )■ e+e^, cu — > tt+tt^ and (f) — > tt+tt^ 
as reported in the Review of Particle Properties (RPP). These are the least experiment 
dependent pieces of informatior0- Actually, as argued in |fT9l (see Subsection 13.3 
therein), we prefer using the RPP information for the products F^ = (1.225±0.071) 10^^ 
and F^ = (2.2 ± 0.4) 10^^ 

We will be even more constraining by supplementing these cu and branching ratios 

'•^Nevertheless, one should keep in mind that these accepted values are highly influenced by the e+e^ tt+tt^ 
scan data samples compared to others. Therefore, this choice could favor the CMD-2 and SND data samples 
within fit procedures; however, as these accepted values are certainly not influenced by none of the BaBar or 
KLOE data samples, the behavior of each of the various ISR data samples becomes a crucial piece of information. 
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by phase information : The so-called Orsay phase concerning the uj decajo and the 
reported phase0of the — t- tt+tt^ amplitude relative to p° — )• tt+tt^. 

• (kk) Use directly data when possible. Indeed, all relevant IB information carried by 
p° — )■ e+e^ and uj vr+vr" can be derived within the BHLS model by the difference 
between the tt+tt^ spectra and the dipion spectrum vr^vr^ in the r decay; more precisely 
using the tt+tt^ spectrum within the tiny energy region 0.76 -r- 0.82 GeV allows to derive 
the relevant IB pieces of information in full consistency with our model. 

As all scan (NSK) e+e" tt+tt- data samples [l25l |26l HH IM 1^ and both KLOE 
data sets (KLOE08 and KLOE 10) stop below 1 GeV, the (p information should be taken 
from somewhere else, namely from the RPR Fortunately, the (j) region is now covered 
by the BaBar data set [|33l[34l . Therefore, as soon as the consistency of the (j) — t- tt+tt^ 
information carried by the BaBar data and by [|24l l48l is established, this part of the 
spectrum could supplement the scan and KLOE data set^ 

Concerning the </> — t- tt+tt^ information, this second strategy will be used by either taking 
the Bill SHI information or the BaBar data points between 1.0 and 1.05 GeV. 

3.3 How to Implement a;/0 — > tt+tt PDG Information? 

The vector meson couplings to tt+tt^ or e+e^ depend on the s-dependent "mixing angles" 
a{s), /3{s) and 7(s). This does not give rise to any ambiguity as long as one deals with spectra; 
however, when using the PDG information for vector meson decays, especially to tt+tt^ or 
e~^e~, one has to specify at which value for s each of the vector meson (model) coupling 
should be evaluated. 

Within the HLS model, there are a priori two legitimate choices for the mass of vector 
mesons; this can either be the Higgs-Kibble (HK) mass which occurs in the Lagrangian after 
symmetry breaking or, especially for the cu and (p mesons, the experimental (accepted) mass as 
given in the RPR Prior to the availability of the BaBar data [|33l , the published e+e~ — t- tt+tt" 
cross section data did not include the (p mass region and, therefore, there was no criterion to 
check the quality of each possible choice in the mass regior0. The choice made in the 
previous studies using the broken HLS model [[TSl [T9l l20l [T3l was the HK mass. 

If one feeds the broken HLS model with the data listed in i-iv (see Subsection 13 . 1 1 above), 
one can make predictions for the pion form factor independently of the measured e+e" — )■ 
Ti~^7i~ data. This procedure is discussed in detail in the next section. Here we anticipate some 
results specific to the mass issue. 

Fig. (d^) displays the prediction for the pion form factor in the region using the HK mass 
to estimate the 07r^7r+ coupling constant with the BaBar data superimposed (not fitted); it is 
clear that the prediction is quite reasonable up to ~ 0.98 GeV as well as above ^ 1.05 GeV. 

'^We will use as input the value 104.7° ± 4.1° found by ll46l . which is consistent with the results recently 
derived [TT] while using an analogous (HLS) model (see Tables VI-IX therein). 

'"^The single existing measurement —34° ± 4° is reported by the SND Collaboration [|48ll . 
Actually, the few BaBar data points between, say, 1.0 and 1.05 GeV carry obviously more information than 
the branching ratio and the "Orsay" phase at the mass. 

'^As the HK mass for the lj meson coincides almost exactly with its accepted RPP value, the problem actually 
arises only for the (p meson. 
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Figure 1: The e+e^ — )• tt+tt^ cross section around the mass together with BaBar data. 
The curve in (a) displays the prediction using the RPP (p decay information computed at the 
Higgs-Kibble mass; the curve in (b) displays the prediction using the PDG decay information 
computed at the experimental mass. In (c) the PDG decay information is replaced by the 
five BaBar data points located between 1. and 1.05 GeV. 



However, it is clearly unacceptable for the mass region in-between. In contrast, using the 
mass as given in the RPP to extract the 07r^7r+ coupling constant from its accepted partial width 
[|24l provides the spectrum shown in Fig. ([B)); this alternate choice is certainly reasonable all 
along the mass region displayed. Therefore, it is motivated to update our former results ni3il by 
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performing the change just emphasize(|3- 

In order to be complete, it is worth mentioning here a fit result obtained by exchanging the 
PDG/SND (j) decay information with the BaBar pion form factor datil^in the range (1.0^1.05) 
GeV. The result, given in Fig. ([1]:), shows that the lineshape of the BaBar pion form factor at 
the (j) mass can be satisfactorily accommodated. 



4 T Predictions of the Pion Form Factor 

4.1 r+PDG Predictions 

As mentioned before, the charged isovector — Ut-ti^-k'^ dipion spectra are not affected 
by7 — p — w — mixing and hence are of much simpler structure. Supplemented by the basic 
p — cu — (p mixing effects which derive from SU (2) and SU (3) flavor breaking, one has a good 
starting point to fix the parameters of the BHLS model to predict the process e+e^ — )■ tt+tt^. 
Specifically, we are using the data including the channels listed in i-iv of Section [3] together 
with RPP information relevant to fix the IB effects affecting the pion form factor. This method 
is named, somewhat abusivel>0 r+PDG. 

Specifically, the IB effects encoded in = Bi{u — t- e^e^) x Br(a; — )■ tt+tt^) and 
Fff, = Br(0 — )■ e+e") x Br(</) — )• tt+tt^), are taken from the RPP (and reminded in Subsection 
13.21) . For the missing phase information we adopt the result from the fit [j46l for the Orsay phase 
of the CO —7- TT+TT^ amplitude and the result from SND [|48ll for the phase of the (p — t- tt+tt^ 
amplitude^. Following the discussion in the preceding Subsection, the model branching ratios 
and phases are computed at the vector boson masses accepted by the RPP. 

The fit returns aprobability of 89.4% with x^/n^o/ = 553.4/596. The fit quality (x^/?^pomts) 
for each of the fitted channels is almost identical to our results in [fT3l (see the last column in 
Table 3 therein). Each of the PDG information listed above contributes by ~ 1 to the total x^- 
It is also worth mentioning that the dipion spectra from [l2Tl l22l l23l are nicely described up to 
a/s = 1 GeV and provide residual distributions indistinguishable from those shown in Figure 
10 of [fT3l. From this fit, one derives the (r+PDG) predictions for the pion form factor which 
can be compared with the various existing e+e^ — t- tt+tt^ data samples. 

The overall view of the comparison is shown in Fig. |2l This clearly indicates that the data 
associated with the channels listed in i-iv, supplemented with a limited PDG information is 
indeed able to provide already a satisfactory picture of the pion form factor as reported by all 
experiments having published e+e^ — )■ vr+vr^ spectra. 

Let us stress that the predicted pion form factor relies on the tt'^tt^ spectra provided by 
ALEPH dill, BeUe [|23l and CLEO [|22l only up to 1.0 GeV. Therefore, the inset in Fig. |I] 
actually shows the extrapolation of the prediction into the spacelike region with the NA7 data 
[|35l superimposed; this clearly indicates that there is no reason to discard the spacelike data 
from our data handling. One should also note that the extrapolation of the prediction above the 

'^We show later on that choosing the cj) HK mass has produced some overestimate of the prediction for and, 
thus, some underestimate of the discrepancy with the BNL measurement |39, 40J. 
'^The "Orsay" phase for both the uj and cf) mesons have no entry in the RPP. 

'''a preUminary version of the present work was presented 1*4911 at the Workshop on Meson Transition Form 
Factors held on May 29-30, 2012 in Krakow, Poland. Some minor differences may occur with the present results 
due to the fact that the SND phase for cf) tt+tt^ was not imposed in the preliminary work. 
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mass is quite reasonable up to ~ 1.2 ^ 1.3 GeV. This may indicate that the influence of high 
mass vector mesons is negligible up to this energy region. 

In order to make more precise statements, let us magnify piece wise the information carried 
by Fig. [21 Thus, Fig. [3]displays the behavior of the various e+e^ — t- tt+tt^ data samples in the 
(0.3 -T- 0.7) GeV energy region. As a general statement the behavior expected from the existing 
data samples looks well predicted by the r+PDG method. A closer inspection allows to infer 
that the CMD-2 and SND data points {i.e. NSK when used together) are well spread onto 
both sides of the predicted curve; this property is also shared by the KLOEIO sample. Even 
if reasonably well described, the KLOE08 and BaBar data samples are lying marginally above 
the r+PDG expectations. 

Fig. IHdisplays the behavior of the various e+e^ — )■ tt+tt^ data samples in the (0.85 ^ 1.2) 
GeV energy region. Here also the predicted curve accounts well for the data behavior. A closer 
inspection tells that the sparse NSK data are well described. The BaBar data are also well 
accounted for all along this energy interval except for the region. As shown by Fig. ^) 
above, this can be cured and one can show that the difference is mostly due to the phase for 
the — )• vr+TT^ amplitude which departs significantl)0 from those provided by SND [48 1. One 
could also note that both KLOE data samples look slightly below the r+PDG expectations in 
this region. 

One may conclude from Figs. [3]and[4lthat our "r+PDG" predictions are in good agreement 
with the data and that a fit using fully these data samples should provide marginal differences 
between all tt+tt^ data sets El]. 

However, the picture becomes quite different in the medium energy region (0.70 -r- 0.85) 
GeV as illustrated by Figure [51 In this region, our r+PDG prediction follows almost perfectly 
expectation from both KLOE08 and KLOEIO data and the detailed lineshape at the u — p 
interference region is strikingly reproduced. Paradoxically, the NSK data are slightly less fa- 
vored - especially around 0.8 GeV - despite their influence on the PDG information used in 
order to account for IB effects in the to — p interference region; however, taking into account ex- 
perimental uncertainties, we already know that a global fit using NSK data is highly successful 

m. 

In contrast, the behavior of the BaBar data looks inconsistent with the r+PDG prediction, 
especially on the low mass side of the interference region. Actually, the observed overestimate 
of the BaBar spectrum affects the whole region from threshold to the to mass but is more im- 
portant in the range (0 . 74 ^ . 78) GeV. At higher energies one observes a reasonable agreement 
with expectations as well as with both KLOE data sets. 

One should also note that the cu mass and (total) width induced by the data for the 7r°7, 177, 
Ti~^Ti~n^ final states are in perfect agreement with all the examined data samples; this indicates 
that the energy calibration around the cu mass is good for all ISR data samples. Fig. ([1]:) has 
already shown that the BaBar energy calibration is also good in the region. 

4.2 Predictions Using the — u Interference Region From Data 

As stated in Subsection 13.11 one can replace the PDG information for to — )■ tt+tt^ and 
p° — )■ e'^e' by any pion form factor spectrum limited to the region (0.76 -r- 0.82) GeV. In 

This issue is examined in detail in Subsection l5.1.4l below. 

Precise information concerning fit qualities is ignored in this Section. 
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particular, this turns out to fit = Bt{ijj — )■ e^e^) x Br(u; tt+tt^) and the Orsay phase as 
it comes out from each of the specified data sets. 

Concerning the tt+tt^ mode, one can also check the BaBar data ll33l [34l versus the 
SND data B8l . This will be discussed below in the context of global fitting (see Subsection 

Using the NSK data or any of both KLOE samples instead of the PDG data does not lead 
to predicted spectra substantially different from their analogue already shown and commented 
upon in the previous Subsection. Some additional parameter constraints have been derived 
from global fits using only the (0.76 -ir 0.82) GeV region from the NSK and KLOE spectra for 
e'^e~ — )■ 7r+7r~; the corresponding results are given in Table[I]and are in reasonable agreement 
with RPP expectations for the branching ratio product and the Orsay phase [46, 47J as well. 

As far as BaBar is concerned, the situation looks different and the most relevant piece of 
information is provided in Figure[6l This proves that the largest difference between BaBar data 
and the other analogous data samples If27l |29l 1311 [321 is the F^^ information inherent to the 
BaBar data. With a (7 ^ 8)a difference with expectations, the BaBar value for F^ is in clear 
conflict with the RPP value, while NSK and both KLOE values are in much better accord. This 
strong disagreement is substantiated by comparing Figures [5] and [6l More precise information 
can be read off Table [T] 



Data Sample 




Local Fit 

Orsay Phase (degrees) 


Prob. (%) 


Input values It24ll46ll 
PDG Refitted value 


L225 ± 0.071 
1.157 ±0.053 


104.7 ±4.1 
108.92 ±2.36 


89.4% 


NSK Hall m 


1.219 ±0.043 


106.71 ±0.25 


90.0% 


KLOE08 Il3ll 


1.076 ±0.041 


110.44 ± 1.24 


87.8% 


KLOE 10 (321 


0.973 ± 0.045 


113.62 ± 1.63 


92.7% 


BaBar [|33l 


1.780 ±0.011 


107.50 ±0.19 


54.4% 



Table 1: Fit results for F^^ = Bt^uj e+e^)Br(a; — > tt+tt^) and for the Orsay phase (in 
degrees) using only the (0.76 ^ 0.82) GeV energy region of each tt+tt" data sample. The fit is 
done following the (local) procedure sketched in Subsection l4.2[ 



4.3 Isospin Breaking Effects in the BHLS Model : Comments 

It follows from the developments just above that the BHLS model fed with a limited num- 
ber of accepted values for some IB pieces of information is indeed able to provide a quite 
satisfactory prediction for the e+e" — )• tt+tt" cross section once the r spectra are considered. 
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This gives support to our breaking model, especially to the s-dependent vector meson mixing 
mechanism. 

The prediction is found in accord with the scan (NSK) data samples and with both KLOE 
data seto- Indeed, the predicted lineshape strikingly follows the central values from both 
KLOE data samples; for the scan data, the prediction based on PDG information is good but 
not as good as for the KLOE data. However, we have shown that changing the PDG requested 
IB information by less than 1 a - as following from a mere comparison of the first and second 
lines in Table [T]- leads to a perfect description of the NSK spectra over the whole available 
energy range. In contrast, the RPP branching fraction product has to be changed by about 
7cr in order to yield a comparable description of the BaBar [13311 data. 

Basically, our approach is a r based prediction of tt+tt^ spectra; it relies on the consis- 
tency of several different physics channels, the r spectra and on a model of isospin symmetry 
breaking. It is thus interesting to examine the consequences of this r based approach on the 
muon g — 2 estimated value. This is what is shown in Figure |7l The first line displays the 
results obtained using the p — )■ e^e^ partial width, the and values extracted from the 
RPP L24J and the two phases already quoted (pure "r+PDG" prediction). The four following 
lines are obtained by replacing the p and uj IB information by the 0.76 ^ 0.82 GeV region of 
the quoted data sets. The line named BNL displays the experimental result [|39l l40l and the 
last line shows the r based estimate from [8 |. The behavior of the BaBar form factor on the 
mass interval [2771,^, compared with expectation explains the difference we get iox g — 2 
compared with [|34l . 

It is clear that all methods used to include IB effects within our r based approach give 
consistent results, all distant from the BNL measurement at the 4cr level. 

5 Global Fits Using the e+e" — > tt+tt" Spectra 

As in our previous analysis [|T3l . we have performed global fits using simultaneously all 
e+e^ annihilation data into the tt^'^, rj'j, Tr+Tr"?!^, K^K^, K^l^ final states, the dipion spectra 
collected in the decay of the r lepton |[2Tll23ll22l and the decay information listed in Subsection 
13.11 We also use the RPP decay properties in the (updated) way emphasized in Subsection 

m 

For what concerns the e+e~ — )■ tt+tt^ data included into the global fit procedure, we have 
performed fits using separately the NSK, KLOE08 and KLOE 10 data samples. Global fits have 
also been performed for the BaBar dipion spectrum restricted to the range of validity of our 
BHLS model. Two options have been considered, using BaBar data up to 1 GeV supplemented 
by —7- TT+TT decay properties from the RPP or using BaBar data up to 1 .05 GeV, thus includ- 
ing the BaBar region and avoiding the need to use RPP information about the 0. In all cases, 
the errors (and the x^) were constructed following the information published/recommended by 
the experimental groups who collected these data. 

We have also performed global fits using combinations of these individual tt+tt^ data sam- 
ples. In this case, the contributions of the NSK and KLOE 10 data to the total were left 
unweighted as their own /n contribution is always of the order 1 in fits using each of them 

^^Some issue with the uncertainties of the KLOE08 data sample will be discussed below. 
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in isolatiorl^. In contrast, in such combinations involving the KLOE08 and/or BaBar data, the 
contribution of each of these to the total was weighted by the ratio fu = n-M / xli (M= 
KLOE08, BaBar) where Xm is the x"^ of the M data set obtained in the best fit using only M as 
e~*^e^ — )■ 7r+7r~ data set; um is the corresponding number of data points. In fits involving the 
spacelike data [|35l [36l. the corresponding weight was also usecj^. 

For definiteness, when relevant, we have used /nsk = fKLOEw = 1^ fKLOEos = 60/90 ~ 
0.67, f BaBar = 270/346 ~ 0.78 and f space = 59/85 ~ 0.69. These weights have been varied 
and it has been found that the sensitivity of the physics results to their precise value is marginal; 
the main virtue of these weights is to provide probabilities not too much ridiculous. 

We have preferred using this method, which turns out to consider each data set as a global 
object, rather than defining local (s-dependent) averages as done by others [6J. This method 
looks better adapted to the global fit method which provides a quality check reflecting the 
behavior of each e+e^ — t- tt+tt^ data set within the global context of a large number of physics 
channels. Indeed, doing local averages would prevent to detect discrepancies originating from 
some given data set only. 

For completeness, it is worth noting that the x^ contributions of the - more than 40 - data 
sets associated with all the other channels were always left unweighted. 

A feature common to all fits using the e+e^ — t- tt+tt^ data sets in isolation or combined 
is that the individual x^ contributions associated with the other channels (7r°7, 777, tt+tt^tt^, 
K+K-, K^K , ...) were only marginally affected by the specific choice of e+e — ^ tt+tt 
data submitted to fit. Their typical values are almost identical to what can be found in the 
last data column of Table 3 in [.13,1 ; more precisely, the x^ value provided by each of these 
channels never varies by more than a few percents. Let us remind that the number of data 
points submitted to fit - beside the e+e^ — t- tt+tt^ data - is ~ 600 when working withiii^the 
configuration B defined in [fT3l (~ 675 within configuration A). 

5.1 Salient Features of the Various e+e~ — > tt+tt" Spectra 

We have performed several tens of fits of the various e"'"e~ — )■ tt+tt" spectra in standalone 
mode and/or combined. It does not look useful to report on each fit in detail. Instead of 
overwhelming the reader with unnecessary information and plots, we have preferred focusing 
on the salient features of their behavior within the global fit context. Beside the fit properties 
of the full spectra (up to 1 GeV, generally), this covers the muon anomalous magnetic moment 
value and the behavior at the u mass, more precisely the value for F^^ = Br(cj — t- e+e^) x 
Br(a; — )• tt+tt"). The first of these topics will be addressed in a separate Section below. 
Concerning the second topic, one reminds that the value for expected from the RPP is 
(1.225 ± 0.071) 10-\ while the r+PDG outcome is rather (1.157 ± 0.053) 10"^ 

On the other hand, we do not enter into much detail concerning the effects of the spacelike 
data, always used weighted in this paper; we limit ourselves to mentioning that they never 

We use the wording "in isolation", "isolatedly", or "standalone" to qualify fits performed using a specific 
e+e^ — > 7r+7r~ data sample. It should be reminded that this tt+tt^ data set is always fitted together with all the 
other data listed in Subsection 13. II unless otherwise stated. 

^"^In the fits referred to in |49|, the spacelike data contributions to the total were left unweighted. 

The difference between configuration B and configuration A is that the former excludes the use of the 3-pion 
data collected around the (p mass. 
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modify the fit qualities in a significant way. 

5.1.1 Standalone Fits of the e+e~ — > tt+tt Spectra 

When using onl>J23]the (unweighted) NSK data, the fit return^ x%skIN^+^- = 128.30/127 
and a global fit probability of 96.3%. We do not show the fit residuals for the scan data which 
are quite similar to Figure 5 in [[TSl . 

Similarly, the KLOEIO data set [32] returns Xkloew/^tt+tt- = 73.68/75 and a global 
fit probability of 87.7%. The fit residual distribution is shown in the top panel of Figure [8l 
the bottom panel in this Figure shows the fractional deviations from the fitting function. Both 
distributions can be considered as reasonably flat. In order to draw the plots, we have used 
as uncertainties the square root of the diagonal elements of the (full) error covariance matrix. 
Therefore, the value for Xkloew '^he flatness of the residuals shown in Figure [8] illustrate 
that the full error covariance matrix is correctly understood. 

Therefore, within the global fit context, each of the NSK and KLOEIO data samples exhibits 
the same outstanding behavior On the other hand, the fit parameters and error covariance 
matrix allows to derive, using obvious notations, F^^^ = (1.205 ± 0.042) 10^^ close to the 
PDG value reported at top of Tabled One also gets F^^oiJio ^ (1,074 ± 0.051) 10"^ a 
~ 2.5 cr difference with the central value for F^^^ . 

Under the same conditions, the (unweighted) KLOE08 data set returns Xkloeos/^tt+tt- = 
96.55 / 60 and a global fit probability of 56%, which is rather low. It is worth noting the remark- 
able flatness of the residual distributions displayed by Figure [9] which does not prevent to yield 
a relatively large value for Xkloeos- The large Xkloeos ^^'^ the flatness of the residual distri- 
bution shown in Figure [9l considered together, might indicate an issue with the non-diagonal 
part of the full error covariance matrix. Anyway, one can conclude that the poor KLOE08 fit 
probability reflects an issue with the KLOE08 error estimate rather than a distorted lineshape. 

Figure[9]should be compared with the similar distribution derived formerly using a primitive 
version of the BHLS model (see Figure 3 in [[T9ll ). The clear improvement substantiates the gain 
provided by the BHLS model in its present form. On the other hand, one should also note that 
the value for = (1.117 ± 0.042) 10"*^ is consistent with F^^^^^^. 

The fit of the (unweighted) BaBar data leads to x\aBarlNn+^- = 343.08/270 = 1.27 
when Umited to 1 GeV (17% probability) and to XBaBar/^n+n- = 340.77/275 = 1.24 when 
going up to 1.05 GeV (corresponding to a 22% probability). For reasons already noted, these 
probabilities are relatively low, reflecting an issue with the rest of the physics involved in the 
global fit, especially the r data. Qualitatively, this result could have been expected from the 
r-i-PDG predictions already discussed. The residual distributions are given in Figure [TOl and 
look reasonable. The residual distribution corresponding to fitting up to 1.05 GeV clearly 
shows a modified behavior in the to region. This has some consequence on the BaBar estimate 
of the muon anomalous moment, as will be emphasized in Subsection 16. 1 1 below. Finally, it 
is interesting to note that the fit outcome provides F^"-^"'' = (1.628 ± 0.012) 10"^ (fit up to 

The agreement of this result with the corresponding one given in Table 3 of ifTSll clearly proves that the 
modification of the (p tt+tt^ information is marginal for the fits not containing the (/> region. The latter is 
covered by the BaBar data set only. 

As emphasized in our previous works liT9ll20]|T3l . the global probabilities are enhanced thanks to the highly 
favorable fit properties of the e+e~ (71^/1])^ data; indeed, /N for these are respectively ~ 0.76 and 0.66. 
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1 GeV) or Fjf"^'^^ = (1.575 ± 0.010) 10"^ (fit up to 1.05 GeV), both being far from F^^^, 
pKLOE08/w fj-Qjj^ expectations lIM . 

5.1.2 Fits Combining the e+e — )• tt+tt Spectra 

In view of the standalone fits just reported, we have done fits of different combinations of 
the scan (NSK) and ISR data samples. In this case, the contributions of the KLOE08 and BaBar 
samples to the (minimized) function are always weighted as already stated. 

• Combining the KLOE08 and KLOEIO data : This returns a consistent picture where 
Xkloeio ^^'^ Xkloeos almost unchanged compared to their standalone values and 
the fit probability reaches 81.6%. This, indeed, confirms that they share the same physics 
content. This is confirmed by the fit resuk i^^io^^os/io ^ ^^^21 ± 0.038) IQ-^, consis- 
tent with both of pKLOEQ?, pKLOElO 

• Combining all the ISR data sets (KLOE08, KLOEIO and BaBar) : This returns xjsR/N^+n- = 
1.34 when including the weights for KLOE08 and BaBar and Xisr/Nt^+tt- = 1-64 when 

the weights are not included. The weighting used does not prevent the probability to re- 
main pooi^(1.8%) reflecting the level of inconsistency of the BaBar and KLOE(08/10) 
data samples already noted. In this case, one gets F^^^ = (1.608 ± 0.010) 10^^, exhibit- 
ing a large distortion towards the BaBar lineshape, despite the weighting. 

• Combining all e+e^ — t- tt+tt^ spectra : Taking the weights into account - which is more 
favorable - one gets X%iobai/^ tt+tt- = 1-35 and a fit probabilit>E3 of 1.3%. Once again, 
the lineshape of the fit function is highly influenced by the BaBar sample in the p — cu 
interference region as shown by = (1.582 ± 0.089) 10^^ 

• Combining the NSK and KLOEIO data : In this case, there is no weight and one gets 
XNSK/^TT+n- = 131.37/127 and x^LOBio/^TT+TT- = 72.90/75 close to the standalone!231 
fit results, i.e. Xkloeio+nsk/^tt+tt- = 1-01, and the remarkable global fit probability 
of 96.9 %, showing that the NSK and KLOEIO data samples are quite consistent with 
each other and with the rest of the BHLS physics as well. In this case, one also yields 

^NSK+KLOEIO ^ (i^igg ^ 03g) 10-6 ^ 

which should supersede the present world average value [|24l because of the full consis- 
tency it exhibits with the largest set of data ever fitted simultaneously and with a splendid 
probability. 

5.1.3 Fits of the e+e — ^ tt+tt Spectra : Concluding Remarks 

The results reported just above have allowed us to show that the NSK data are in fair agree- 
ment with the physics represented by the tt°7, 777, tt"''tt"tt°, K^K~, K^K^ annihilation chan- 
nels, the dipion spectra collected in the decay of the r lepton and some more decay information 
listed in Subsection 13. II This is not a really new result as this conclusion was already reached 
in our [|T3l . The single difference with [fT3l is the new dealing with the — t- tt+tt^ information 
(see Sub section [33]) imposed by the BaBar data. 
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The new information is that the KLOEIO data sample behaves likewise and, more impor- 
tantly, that the NSK and KLOEIO data sets are consistent with each other as well as with the 
rest of the physics considered within the BHLS model and the global fit context. 

We have performed global fits which have shown that the data from KLOE08 and BaBar 
have some difficulty to accommodate the global fit context. For what concerns KLOE08, the 
problem appears to be related with (underestimated?) systematic errors or, possibly, with their 
correlations. In the case of BaBar data, the problem looks more serious, as it deals with the 
form factor lineshape itself in the p — to interference region; this issue manifests itself in the 
value for F^, much larger than expected. 

For this purpose, it is interesting to see the behavior of the most relevant fits of the pion 
form factor in the p — to resonance region. These are displayed in Figure \TT\ The best fit 
obtained using only the NSK data is shown in the top left panel and is clearly quite satisfactory. 
The top right panel exhibits the case when the NSK data are complemented with the KLOEIO 
sample; this fit is also quite successful. 

Bottom left panel in Figure [TT] shows the behavior of the best fit function when only con- 
sidering the ISR data (BaBar, KLOE08, KLOEIO) and bottom right panel when taking into 
account all existing scan and ISR data. Both are clearly less satisfactory reflecting mostly the 
tension between BaBar and the KLOE data sets. 

The question is now whether the KLOE08 and BaBar data samples can nevertheless help 
to improve some physics information of important concern as the muon g — 2. This will be 
discussed in the forthcoming Section. Anyway, the KLOEIO data sample allows to confirm the 
results already derived using the NSK data and, even, helps in getting improved results. 

5.1.4 The (j) Region in the e+e tt+tt Spectrum 

Up to now, the pieces of information used in our fits/predictions are the RPP value for the 
product F,^ = Br(0 e+e~)Br(0 tt+tt") = (2.2 ±0.4) 10"^ and the "Orsay" phase for the 
(p —7- TT+TT^ amplitude provided by SND [48], namely —34.0 ± 4.0 degrees. In order to avoid 
over interpreting the SND phase as an Orsay phase (i.e. identified with the phase of the product 
F^^{s) 5'(/,7r7r(s) in Eq. (fTTI) at s = m|), we have found worth revisiting this assumption. Using 
the NSK and KLOEIO data sets as reference vr+vr^ data sample, we have performed several 
global fits within the BHLS framework. 

For definiteness, and for an easy comparison, the first line in Table |2] displays the values 
extracted from [|24l l48l and used in all fits referred to up to now. The second line in this 
Table displays the corresponding values reconstructed using the final parameter values from 
the appropriate global fit. The goodness of fit (94.5% probability) is reflected by the values 
shown in the last two data columns. The reconstructed value for is in very good accord with 
the (input) RPP value, while the phase is shifted toward slightly more negative values than the 
SND datum. 

The fit results summarized in the third line are obtained by withdrawing the SND phase 
||48l from the minimization procedure. The fit quality is almost unchanged (95.3% probability) 
and the expected is still well reproduced; on the other hand, the phase moves by about 3a 
but remains significantly negative. 

In both cases, the pion form factor lineshape exhibits the behavior shown in Figure ^Bp), 
e.g. a tiny peak at the mass, resembling what is shown in Figure 4 of [i48il . 
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The next step has been to remove the reference values for and for the phase and 
supplement the NSK and KLOEIO data with the region spectrum from BaBar data (y/s E 
[1.00, 1.05] GeV). The important results are shown in the fourth line of Tabled The fit proba- 
bility remains good (91.3%), but the of the NSK+KLOEIO data is degraded by ~ 10 units, 
which may indicate some tension between these data sets and the coherent background beneath 
the peak. The fit value for increases by 50% but remains perfectly consistent with expec- 
tations [l24l. In contrast, the Orsay phase changes dramatically (from —[30 ^ 50]° — )■ 150°). 
The lineshape exhibited by the BaBar pion form factor is almost identical to those shown in 
Figure ([lb). 



Information Used 


F^IO^ 


Orsay Phase 


[Xn+TT~ /-^points] <1.0 GeV 


[Xtt+tt- / -^points] <j) 


Input 11241 1481 


2.2 ±0.4 


-34.0 ±4.0 






Imposing Orsay phase [|48]| 


2.34 ± 0.42 


-41.45 ± 1.89 


204.3/202 


2.25/2 


Releasing Orsay phase 


2.23 ±0.41 


-48.23 ± 1.88 


202.9/202 


0.04/1 


BaBar Region Spectrum [[33ll 


3.31 ±0.99 


156.92 ±2.76 


215.0/202 


5.2/5 


BaBar data |l33l fit 
Imposing BaBar parameters 


3.04 ±0.60 

3.05 ±0.50 


153.41 ±0.98 
153.37 ±0.85 


(338.1/270) 
208.0/202 


(2.69/5) 
0.05/2 



Table 2: Fit results for = Br(0 — )• e'^e )Br(0 — )• tt+tt ) and for the "Orsay" phase (in 
degrees) for the amplitude. See Subsection 15 . 1 .41 for comments. 



As stated abovJEl, the [1.00, 1.05] GeV region of the BaBar spectrum carries information 
on the signal but also on the underlying (p coherent) background. In order to subtanstiate 
its effect, we have finally performed a fit using as input the and phase values extracted 
from fitting the BaBar data up to 1.05 GeV (see Subsection 15. 1.11) : the corresponding values 
are given in the fifth data line of Table [2] and, when imposed to the fit procedure, these numbers 
lead to the results shown in the last line. The fit probability improves to 94%. 

As a summary. Table [2] shows that some ambiguity occurs about the pion form factor be- 
havior in the region. It becomes interesting to have new data covering this region, in order to 
decide which among Figures ([IJ5) and (dt) reflects the right behavior. It is also interesting to 
examine the consequences of this ambiguity on g — 2 estimates. 

Without going into details (addressed in the next Section), we can give the values for Aa^ = 
a^^P — a*^ corresponding to the various cases shown in Table [2l In units of 10^^°, the values 
for Aa^ are respectively 39.91 ± 5.21 (line # 2), 39.35 ± 5.23 (line # 3), 37.18 ± 5.16 (line # 
4) and 39.32 ± 5.20 (line # 6). The most likely effect this implies on Aa^ is a possible shift 
by 0.59 X 10-^° (lines # 2, 3, 6). A maximum possible shift is 2.73 x 10"^° (lines # 2 and 
4). Instead, the uncertainty is practically unchanged. Taking into account the uncertainty on 
^exp [|39| [40|| ^ the significance implied by this shift is quite marginal (from 4. 880" to 4.57(j). 
Nevertheless, this shift plays as a systematic effect to be accounted in the final value which will 
be proposed for a*f . 
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6 Hadronic Contribution to the Muon g — 2 Estimates 



In order to estimate the hadronic contribution to the muon anomalous magnetic moment 
a^, the method followed in the present work is identical to the one used in [|T3| . In this Section, 
we examine the global fit solutions provided by the data samples listed in Subsection 13. 11 the 
various cases correspond to varying the e+e^ — )■ tt+tt^ data sample combination submitted to 
the global fit procedure. 

In our previous analysis [|T3l . we encountered some tension in the e+e" annihilation data in 
the region between the e+e~ 7r+7r~7r° channel and the e+e^ K^K^, K^K^ channels. 
Including all data in the global fit defined what we called "configuration A". Alternatively, 
we excluded the 3-pion data in the vicinity of the (p mass from the fit. In this case, the cor- 
responding collection of data included in the fit was called "configuration B", which yields a 
better fit probability. As [fT3l , the present analysis privileges configuration B as a basis. As the 
"observed" tension could not be abnormal, we also consider fits based on configuration A for 
completeness. 

6.1 a^{'KTT) Contribution from the Reference mTrvr Region 

In order to emphasize what is going on, it is worth examining the contribution to a^(7r7r) 
provided by the invariant mass interval A = [0.630, 0.958] GeV. As for this reference region, 
all experimental groups have published the numerical integration leading to their estimate of 
a^iy-K-K, A) (using Eq. ([T]) with the experimental spectra), one can compare our fit outcome with 
these. When combining data sets within the fit procedure, we compare to the usual weighted 
average of the a^(7r7r, A) values provided by each experimental group. 

The results for a^iixu:, A) are collected in Fig. [I2land will be discussed from top to bottom. 
The first point has no experimental partner, as it comes out from the r+PDG fit described in 
Section IH The dashed-dotted line is drawn through the central value of this prediction, the 
dotted lines show the corresponding ilc band. 

It should then be noted that most experimental estimates happen to be inside the la band 
of the r+PDG prediction. KLOEIO is at the border, while BaBar is rather larger by about 
2cr. One should also note that KLOE08 as well as the combination CMD2+SND+KLOE10 
(our preferred combination because of its global fit properties) provide central values which 
coincide almost exactly with the r+PDG expectation. This looks to us a noticeable property 
taking into account the long-standing discussion about the r — e+e^ issue. Our results tend to 
prove that this issue vanishes once the breaking mechanism is appropriate. 

The next point shows the global fit estimate (square symbol) while using solely the NSK 
data (especially those from [|27ll28ll29l , but also those from [|25ll26l ). The present fit value is 
found smaller than those published in [fT3l by about 1 x 10~^°. This shows that the effect of 
our updating the (j) data decay information extends over the whole spectrum. The uncertainty is 
unchanged and amounts to a 40% gain compared with the experimental estimate. This gain is 
a pure effect of the global fit procedure where the channels others than tt+tt^ allow to improve 
also the tt+tt^ contribution importantly because of the underlying VMD physics correlations. 

The following point is derived using solely the KLOE08 data sample. As with fitting solely 
the NSK data, the fit outcome also differs by 1 x 10^^° from the experimental central value. 
The fit result associated with the KLOEIO data sample is also found in good agreement with 
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the experimental expectation and one should note that both differ by 1.5 x 10^ in the opposite 
direction. In contrast, the fit result for the mixed KLOE10+(weighted) KLOE08 data sample is 
found to agree very well with expectation (a 0.6 x 10^^° difference only). 

The data point associated with the CMD2+SND+KLOE10 combined data sample is found 
in much better accord with the average of the experimental data; in this case, the error is also 
improved by about a factor 2 compared with the experimental value, but, interestingly, first, the 
central reconstructed value does not exhibit any significant systematic shift (~ 0.1 x 10"^") 
compared to the experimental expectation, second, the central reconstructed value is in perfect 
accord with the r+PDG prediction, as already noted. 

The next pair of points shows the case for the BaBar data. While using the whole spectrum 
up to 1 GeV, the fitting function provides an estimate for a^^mr, A) in good agreement with 
the expected experimentally integrated value. However, extending the fitting region up to 1.05 
GeV, increases the difference between the fit outcome and the experimental expectation from 
1.5 X 10~^° to 4.4 X 10^^°, which illustrates some issue when fitting also the BaBar region. 

The last two points (ISR and scan + ISR data, resp.) correspond to fitting the whole spectra 
up to 1 GeV while weighting the KLOE08 and BaBar contributions to the x^- The results look 
in good correspondence with the weighted averages of the direct integration results, and one 
remains within the la band of the r+PDG expectation. 

The contribution to a^(7r7r) from the invariant mass interval A = [0.630, 0.958] GeV 
and the global fit properties (see Section [5]) allow to draw a few important conclusions : 

• (1) There is no mismatch between the r+PDG expectations and the estimates derived 
from the fitting functions or those derived from the numerical integration of the measured 
e"'"e^ — 7- TT+TT^ cross sections. To be more precise, when some departure is observed, 
it is always closely associated with poor global fit qualities of the corresponding data 
sample. 

• (11) One can take as reference e+e" tt+tt" data the NSK (CMD-2 & SND) and 
KLOEIO data samples. Separately and together, they are found in perfect accord with 
all annihilation data considered (tt^'J, rj'j, tt+tt^tt^, K^K^ , K^l^ final states) as well 
as with the published dipion spectra in the decay of the r lepton and the few additional 
decay information introduced within the global fit BHLS framework. The uncertainty is 
improved by a factor of about 2. 

• (111) Comparing the central value for a^{mi, A) derived from the fitting function - using 
the NSK and KLOEIO data samples within the global framework - with its direct esti- 
mate indicates that the estimate derived from the fit is almost unbiased (~ 0.1 x 10^^°). 
In contrast, the fitting functions derived from the fits to each of the NSK and KLOEIO 
data separately may exhibit a small bias (a shift of ~ 1 x 10^^°). 

• (Iv) One may remark that the global fit using all existing e+e^ — )► tt+tt^ scan & ISR 
data also provides reasonable values for a^{7ni,A) without a significant bias and with 
a much improved uncertainty. This reduced uncertainty may look too optimistic, taking 
into account the global fits properties of the corresponding sample combinations. 
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6.2 The Full Muon Anomalous Magnetic Moment a 



Here we are interested in the theoretical prediction of the muon anomalous magnetic mo- 
ment and its comparison with the experimental result [39^, '40'|. Our BHLS based global fit 
results provide the dominant contribution to the hadronic vacuum polarization contribution. 

Several global fits have been performed using different groups of e+e^ — t- tt+tt^ data sam- 
ples. The resulting fitting functions and the parameter error covariance matrices have allowed 
to derive the leading order hadronic contribution (LO-HVP) to up to 1 .05 GeV in various 
cases. The full HVP is obtained by combining our present results with the other hadronic con- 
tributions listed in Table 9 of our previous work [|T3l . The various evaluations of the muon 
g — 2 are calculated by summing up a'^^^^^^ with the Higher order HVP, the QED and the 
electroweak (EW) contributions and the so-called light-by-light (LBL) contribution; for this 
purpose, the values given in Table 10 of |fT3l have been updated (see Table [3] below). 

In Figure [131 we display results for the deviation Aa^ = a^'^P — a}"^^ between experiment 
and theoretical predictions in units of 10^^°, for various fit results which only differ by the 
indicated e+e^ — )■ tt+tt^ data sample considered within the fit procedure. 

Specifically, the fit results included in Fig. [T3]are all based on the "configuration B" data set 
introduced in [|T3il . Furthermore, always included in the fits are all channels listed in Subsection 
13. 1[ especially the r data. 

At top of the Figure, we provide the r-i-PDG entry {i.e., without using any e^e~ — )■ tt+tt^ 
data set in the fit procedure). The twin points drawn correspond to Aa^ = a^^^ — a^^ in 
units of 10^°; the lower point is obtained using the spacelike data |I351[36]|, the upper one by 
excluding them from the fit. The number following the point is the estimate for Aa^ from fits 
excluding the spacelike data. The significance of this difference is written [x/y a], xa being 
the significance obtained when excluding the spacelike data, ya when they are included into 
the fit procedure. The last column in this Figure provides the /upoints for the e+e^ — )■ tt+tt" 
data; this number is given only for (combination of) samples used without any weight. 

The first remark which can be drawn is that including the spacelike data does not change 
noticeably the value for Aa^ nor its statistical significance. The second important point is that 
all the features noted while discussing a^(7r7r, [0.630, 0.958] GeV) in the previous Subsection 
survive when dealing with the whole provided spectra (see (l)-(lv) just above). 

It was noted in lfT3ll that the use of the ISR data should not dramatically improve the un- 
certainty of the theoretical estimate : going from scan data only (the NSK entry) to scan-i-ISR 
data (the "Global Fit" entry) improves the uncertainty by only 0.3 x 10^^°; this is the major 
gain of working in a global fit context [rT3l|. In contrast, the central value exhibits slightly more 
sensitivity, as can be seen by comparing the various entries. 

Anyway, the present analysis allows to conclude that the discrepancy between the measure- 
ment for 1391 |40ll and the Standard Model prediction is certainly in the range of 4.5(7 or 
slightly higher. Taking into account the already reported fit qualities, the most favored value 
for ajo-^^^ is : 

(j^LO-HVP ^ 23 ± 4.51) 10"^°. (13) 

This value is derived by introducing the e+e~ NSK and KLOEIO data samples within our 
global fit procedure. 

Finally, we have redone several of our fits within the configuration A. These confirm the 
results derived within the configuration B which are the body of the present paper. One also 
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confirms the 'tension' reported in |fT3l| between the 3-pion spectra in the (j) region and the data 
collected for both e+e^ — )• KK cross sections. Whether this is an experimental issue, or if this 
indicates that our breaking model should be extendecj^ is an open issue. Nevertheless, some 
physics results derived within the configuration A will be presented just below. 

6.3 Comparison with Other Estimates of 

Fig. [l4]reports on some recent estimates of the muon anomalous magnetic moment to- 
gether with the BNL average value [|39ll40l . In this Figure, our (favored) results - derived using 
only the CMD-2, SND and KLOEIO data samples - are given under the entry tags BHLS::A 
(for configuration A) and BHLS::B (for configuration B). Some evaluations proposed by other 
authors are also shown; they are extracted from L6J (DHMZIO), LZJ (JSll), L8J (DHea09) and 
[111 (HLMNTll). 





Values 
scan only 


(incl. r) 

scan © ISR 


Values ( 
scan © ISR 


;excl. r) 

scan only 


LO hadronic 


685.66 ±4.54 


688.60 ±4.24 


687.53 ± 4.34 


684.25 ±5.15 


HO hadronic 


-9.97 ±0.09 


LBL 


10.5 ±2.6 


QED 


11 658 471.8851 ±0.0036 


EW 


15.40 ± O.lOhad ± 0.03Higgs,top,3-loop 


Total Theor. 


11 659 173.48 ±5.23 


11 659 176.42 ±4.89 


11 659 175.35 ±5.06 


11 659 171.89 ±5.77 


Exper. Aver. 


11 659 208.9 ±6.3 


Aa^ 


35.43 ±8.14 


32.48 ±8.03 


33.56 ±8.08 


36.89 ±8.54 


Significance {na) 


4.33(T 


4.05(T 


4.16(T 


4.32a 



Table 3: The various contributions to lO^^O;^. Aa^ — a'^fP — a^'* is given in units of 10^^*^. These resuhs have 
been derived using all existing annihilation data samples (configuration A). The naming "scan" and "ISR" in the 
subtitles refers only to the tt+vt^ channel. For the measured value a^^'', we have adopted the value reported in 
the RPP which uses the updated value for A = Hfj,/fip recommended by the COD ATA group [,501 . 



The analyses reported in [|9|| and |l7l provide g — 2 evaluations based on all existing data 
samples. In particular, all tt+tt^ data samples are used (ISR and scan) and all other existing 
annihilation data samples, especially the 3-pion data samples. These evaluations of g — 2 

^^One might have to revisit the possibility that nonet symmetry should also be broken in the vector meson sector 
as it is already done in the pseudoscalar sector. 
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may also include or not the r data. Therefore, it is also instructive to give the results derived 
from BHLS under the same conditions, i.e. working in the so-called configuration A [13J and 
using all available ISR data samples beside the scan (NSK) data, whatever are the global fit 
probabilities. The corresponding results are reported at the "Full BHLS" entries in Fig. [14] 
More precise information is also provided in Tabled However, one should note that, while the 
fit probabilities are good when using only the scan (NSK) data (90.6 % or 72.7 % respectively), 
they are below the 10 

-4 leveiSH when all the ISR data samples are included. This reflects the 
statistical inconsistency between both KLOE data samples on the one hand and the BaBar data 
sample on the other hand. 

We have explored this issue in cases including the r data samples. Without going into much 
details, it looks worth to comment about fit qualities. The just reported probabilities bounded by 
10~^ are worrying, taking into account that the problematic data samples have been weighted. 
For this purpose, we have first performed a fit discarding the BaBar data sample, to check the 
level of statistical inconsistency of NSK, KLOEIO and (weighted) KLOE08; the obtained fit 
probability is 61.6 % and the numerical results are almost identical to those already derived 
using only NSK and KLOEIO. In particular, the scan data yield x^{.NSK)/n = 136/127 = 
1.07 and F^^ = (1.11 ± 0.03) 10^^, i.e. a limited distortion of the scan data description. 
The same exercise performed using NSK and (weighted) BaBar yield 12.4% probability but 
X^{NSK)/n = 171/127 = 1.35 and = (1.52 ± 0.04) 10^^ i.e. a severely degraded 
description of the scan data. 

It should be noted that, when using all tttt data samples, the value returned by the BHLS 
model is close to JSl 1 An interesting point here is that [7 J uses a mixing procedure (7 — 
p — to) parent to ours {p — to — (p), while the IB effects are introduced differently in order to 
get the DHMZIO result Comparing the Full BHLS values derived when using or not the r 
spectra indicates that BHLS is much less sensitive to the use of the r spectra that the method 
from [6] (compare the DHMZIO entries). Finally, comparing DHea09 [8] and BHLS when 
no ISR data are considered, one observes a good consistency; this seems to indicate that the 
problem lays in the amount of isospin breaking compared with r. 

Comparing now aU the BHLS results (BHLS::A,BHLS::B and FuU BHLS) indicates a re- 
markable stability for Aa^ = a^^^ — a|f under various fit conditions. One should also note that 
the r+PDG prediction is always in quite nice agreement with using any kind of configuration 
for TTTT data samples (BHLS::A,BHLS::B and Full BHLS), including those discarding the r 
spectra and the ISR data (see also Table [3]). 

Nevertheless, one should stress that, because of the poor fit quality yielded when using 
the KLOE08 and BaBar data samples beside the NSK and KLOEIO ones, we do not consider 
the results shown in Table [3] and in Fig. [141 under the entry tag "Full BHLS" as reliable. As 
a general statement, we prefer relying on estimates using only the NSK and KLOEIO data 
samples for which the probabilities are good (~ 90% -r- 98% for BHLS::B, ~ 56% 73% for 
BHLS::A). Numerically, the BHLS::B and BHLS::A estimates happen to be always close to 
each other. Nevertheless, we prefer favoring the BHLS::B result over BHLS::A for reasons 
already explained. 

The statistical significance for Aa^ = a^^^ — a^^ is displayed on the right-hand side of the 
Figure for each of the reported analyses. The updating of the —t- tttt treatment within our 
minimization code turns out to increase the significance for Aa^ by about 0.5cr without any 
significant change in the fit qualities. The present study, however, has shown that this update is 
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highly favored by the information carried by the BaBar spectrum. 

On the other hand, one should note that using all ISR data samples only leads to a ~ 4% 
improvement for the uncertainty on a^^ compared with discarding all of them. This is actually, 
as reminded in the Introduction, a specific consequence of the global fit procedure. 

Before closing this Section, it is worth remarking that the central values for a^^^^^^ 
derived using scan data only or scan © ISR are very close to each other. This is due to a 
numerical conspiracy of the BaBar and KLOE data samples : The former data tend to increase 
a^^~^^^ and the latter to decrease it, by a quite comparable amount. Taking into account the 
above reported fit probabilities, the increase of a^'^~^^^ is not justified on statistical grounds; 
in contrast the decrease of ajf^^^^ produced by KLOE 10 (and also KLOE08) is supported 
by quite good probabilities of simultaneous fits of scan and KLOE data. 



6.4 A Few properties of the Favored (NSK+KLOEIO) Estimate For 

It follows from the present study that a consistent data sample can be defined which allows 
to yield a reliable estimation of . This sample contains all the collected data samples cov- 
ering the 7r°7, 777, K^K~ and K^Ks annihilation channels and the dipion r spectra up to 1 
GeV. It also contains the tt^tt+tt^ annihilation channel. For this last process, we prefer working 
within the so-called configuration B by excluding the data sets collected around the mass. 
However, as clear from Figure O taking them into account (configuration A) does not lead to 
substantial differences. Concerning the fundamental tt+tt^ annihilation channel, the statistical 
properties exhibited by the fits lead us to rely on using only the scan (NSK) and the KLOE 10 
data samples. 

This leads us to conclude that the most favored value for Aa^ = aff^ — aff* is 39.91 ± 5.21 
(in units of 10~^°). S ub section 16 . 1 1 tends to indicate that this value should be almost unbiased. 

However, the ambiguity about the choice of the Orsay phase discussed in Subsection [5X41 
may also indicate a possible systematic shift, likely of 0.59 x 10^^°, but certainly limited to 
2.73 X 10^^°. This upper bound supposes that the BaBar spectrum behavior in the region 
should be trustecj^. Only a devoted measurement of the pion form factor in the region may 
allow to conclude about this possible systematic shift. 

Looking at the differences between BHLS::A and BHLS::B, one also observes that the most 
significant difference between both solutions is concentrated - as could be expected - in the 
contribution of the 3-pion data around the peak. This may lessen Aa^ by 1.26 x 10^^°. 

On the other hand, it is worth noting the behavior of the predicted pion form factor in the 
region rriT,.,, < 0.5 GeV down to the negative s region. Figure [15] shows that the pion form 
factor derived from fits describes quite well the (highly constraining) spacelike data down to 
about —0.15 GeV^. Whether spacelike data are included within the data sample submitted 
to fit (full curve) or not (dashed curve) does not make any difference (both curves coincide 
within the thickness of the lines). This is a noticeable property as, following Section 6 in [[T9l . 
the NA7 data should be rescaled. The rescaling factor 1 — A is such that A depends on the 
spacelike data, on their scale uncertainties (0.8 % for NA7 [35]) and on the fitting function. 
Therefore, by ascertaining the threshold behavior of the pion form factor. Figure [T5]proves that 
the contribution of the tttt threshold region to — 2 predicted using the NSK and KLOEIO data 

^'One may remark that the pion form factor fitting function used in |34| does not include a term. 
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is reliable and confirmed by the existing spacelike data. Indeed, these are well accounted for 
without any regard to their being included within the data sample submitted to fit (see Figure 

est)). 

Anyway, the statistical significance of Aa^ for the preferred solution is 4.88cr. On the other 
hand. Figure [14] has shown that using or not the r spectra does not significantly modify the 
significance for Aa^ (5. Oct — )■ A.9a). In contrast. Figure [14] also shows that using the KLOEIO 
data sample increases the significance by about 0.5a. 

In view of all these considerations, especially the issues encountered with the region 
physics information, the most reliable estimate we can propose is : 



Ah 



Aa, 



(11 659 169.55 + 



+1.26 
-0.59 



± 5.21^;,) 10 



-10 



^exp 



^th 



(39.35 + 



+0.59 
-1.26 



±5.23,ft ±6.3 



Jexp, 



10 



-10 



(14) 



a 



LO-HVP 



(681.82 



f 1.26 
-0.59 



±4.5U) 10 



-10 



The quoted is an estimate of the possible uncertainties affecting the (f) region and 

emphasized in Subsection 15.1.41 and also just above; the central value is derived from a fit 
within configuration B and avoiding the use the SND phase constraint, [^a^'*]^ should not be 
added in quadrature to the theoretical error but linearly to the central value. Accounting for this 
possible shift slightly lessens the significance for Aa^ from A.Sla — )■ 4.67a. 

Therefore, Figure [14] illustrates that the global fit approach allows for a gain in the signif- 
icance for Aa^ of about la or more compared to the traditional methods based on numerical 
integration of the experimental cross sections [[8l[6ll71l9l. It has also been shown that the central 
value for a'^^ coincides almost exactly with its r+PDG prediction. Moreover, our approach has 
allowed us to define the largest set of data samples which exhibit full consistency with each 
other. 

The results involving the fit of the scan © (all/selected) ISR data, do not lead to a noticeably 
smaller uncertainty for Aa^ compared to using only scan data, as already inferred in (13] . In- 
deed, Figure[T4]indicates that in changing the fit conditions [scan ©r] — )• [scan ©r© KLOEIO], 
the uncertainty marginally improves (5.28 5.21), while the central value moves by 3 units 
(36.88 — )• 39.91), increasing the significance (4.5a — )■ 4.9a). 



7 Conclusion and Perspectives 

The e~^e~ — )► n~^7r~ annihilation channel has been widely studied and one is faced with 
several data sets collected by different groups under various conditions. Beside the scan ex- 
periments performed by the CMD-2 and SND collaborations which have produced valuable 
data samples [I26ll27ll28l l29l, experiments using the Initial State Radiation (ISR) method have 
also been performed by the KLOE and BaBar collaborations and have produced (much) higher 
statistics data samples [|31i |32[ [33l . The noticeable gain in statistics - for what concerns the 
vr+vr" intermediate state - is partly balanced by the issues raised by the dominance of system- 
atic uncertainties. Indeed, comparing the properties of these data samples reveals inconsisten- 
cies which leads to somewhat contradictory conclusions concerning the predicted value for the 
muon anomalous moment. 
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It thus becomes a relevant challenge to find a tool able to examine critically the properties of 
the various available samples and substantiate their differences. The present study has proved 
that the existing data for : 

• The e^e^ annihilation to the 7r°7, 777, tt+tt^tt'^, K+K-, K^K final states, 

• The dipion spectra in the — )■ tx^t^'^i' decay, 

• Some radiative decays (namely, tx^ /rf/rj' — )■ 77, — )■ 7r^7, K* — )■ K7, rj' — )• oj'-^, 
(j) — 7- ?7'7) as given in the RPP, 

supplemented with some limited isospin breaking (IB) information, provide a benchmark able 
to reconstruct with a noticeable precision the pion form factor measured in the e+e^ — )• tt+tt^ 
annihilation. As could have been expected, the requested IB pieces should cover the — t- 
TT+TT^ decays, but also the p° — )■ e+e^ decay, rarely stressed explicitly [7]. 

The tool for this prediction is the broken HLS model (BHLS) defined and studied in [T3l. 
In this last reference, the scan data for the e+e~ — )■ tt+tt" annihilation channel were used. 
However, the present study has proved that replacing these data by only the decay data 

just referred to, allows BHLS to pin down all the parameters of the model and provide strik- 
ingly precise predictions for e+e^ — )■ tt+tt^. These have been named "r+PDG" predictions of 
the pion form factor. As these predictions are clearly a new way to formulate r based predic- 
tions for the HVP, BHLS provides also a new tool to explore the long reported r versus e+e^ 
discrepancy between the HVP evaluations. 

Then, the present study leads to the following conclusions : 

• 1/ There is no mismatch between r based predictions and direct e+e~ evaluation of the 
HVP if one relies on BHLS and, especially, on its isospin breaking scheme. This is essen- 
tially implemented in two steps : direct breaking at the HLS Lagrangian level followed by 
vector meson mixing. The vector meson mixing, which generates physical vector fields, 
is unavoidable because, at first order in breaking parameters, the ideal vector fields are 
no longer mass eigenstates. This simply confirms the study in [fT3l . 

• 2/ The r-i-PDG predictions compare astonishingly to the KLOE data samples all along 
the spectra. Because of their influence on the PDG information, it is somewhat para- 
doxical to observe a difference between the NSK data and the r-i-PDG predictions, more 
marked than with both KLOE data samples. However, a la modification of the width for 
oj — ^ TT+TT^ allows to recover a good agreement between the NSK data and the r-i-PDG 
predictions. In contrast, the predicted pion form factor reveals an important disagreement 
with BaBar data, especially in the [0.74, 0.78] GeV region. 

• 3/ By varying the fit conditions, it has been found that the KLOE 10 and the scan (NSK) 
data samples together lead to global fits with outstanding statistical properties. Even if 
the KLOE08 data sample has properties similar to KLOEIO (see Figures |9] and [8]), it 
looks premature to include it within the global fit procedure^. 

^"An extraction of the pion form factor from the KLOE 2008 data is presently under way, which rehes on using 
the 7r7r7/^/i7 event ratio 1,5 U . As this implies a deaUng with uncertainties slightly different from 1,31] . this might 
change the picture. 
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From our analysis, one also gets : 

= Bi{uj e+e") x Br(w ^ Tr+Tr") = (1.166 ± 0.036) 10"^ , 

close to the presently accepted 1241 value and twice more precise, with a corresponding Orsay 
phase of 104.73° ± 0.63°. This value is in accord with the values derived from separate fits 
to scan, KLOEIO and KLOE08 data samples, reflecting that the lineshape of their pion form 
factors are consistent with each other. In contrast, the value derived from (global) fit to the 
BaBar data sample is different by about (7 ^ 8)cr. 

Concerning the meson region, our analysis of BaBar data provides Br(0 — )■ e^e^) x 
Br(0 —7- TT+TT^) = (3.31 ± 0.99) 10^^ in good agreement with the accepted value ll24l . The 
corresponding Orsay phase, however, seems to disagree with expectations ll48l : our fits tend to 
favor -48.23° ± 1.88°, closer to the SND phase. 

Using the KLOEIO and the scan data samples leads to the most probable value for the muon 
anomalous moment : 

af = (11 659 169.55 + 

which exhibits a significance for Aa^ = a^^^ — at a (4.7 -i- 4. 9) a level, significantly larger 
than the results fully derived by direct numerical integration of the experimental cross sections. 
Some additional remarks are worth being made concerning this result and our approach : 

• 4/ The cross sections which should be integrated within the BHLS model to evaluate 

can be considered as an "effective field theory induced interpolation" between data 
points. By using a relatively fine energy binning, the interpolation uncertainties are cer- 
tainly minimized in all regions, and for all cross sections, which exhibit sharp energy 
variations. This is certainly an advantage over the standard method which should rely 
on trapezoidal estimation between relatively distant measured data points, possibly im- 
proved by taking somewhat into account the local curvatures. 

• 5/ One has certainly noted that the value for a^^ we privilege is smaller than all estimates 
involving all ISR data. As explained above, this is because the effects of the KLOE 
sample is not balanced by the BaBar data. However, as noted several times above, the 
value for F^^ is a criterion which leads us to conclude that NSK and KLOEIO only have 
homogeneous properties which justify to privilege a global simultaneous fit of these. 

• 6/ The poor description of the BaBar data within the global fit framework does not mean 
that it cannot yield a good description in a really standalone fit. However, the global fit 
method, which accounts for the underlying physics correlations, performs quite success- 
fully with other data samples covering the same physics and the same energy region. As, 
it is unlikely that BHLS may have some inherent reason to exhibit some tropism towards 
the scan or KLOE data, we indeed consider the ajf value we propose as the most moti- 
vated one on statistical grounds. Of course, new ISR data one can expect from Belle are 
clearly valuable. 



+ 1.26 
-0.59 



± 5.21,.) 10 



-10 
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Nevertheless, the results derived using all tt+tt^ ISR data samples have been examined 
(see Figure [T3l) and shown to agree with the favored estimate just quoted. Because of the 
properties expected from the (BHLS) global framework recalled in the Introduction, using 
only the KLOEIO and NSK data samples provides already as precise results as those derived 
using all available ISR data samples and the traditional evaluation method of a^^~^^^ . 

Our work tells that the significance for Aa^ starts to be close to the ha level. How close to 
this value it could be, should be confirmed by more precise annihilation data in the [0.95 , 1.05] 
GeV mass region, especially in the tt+tt^, tt+tt^tt^ final states. Finally, new measurements 
for the experimental value of the muon g — 2 are planned [|52l[53ll : they are important so as to 
confirm the central value for a^^^^, and also to lessen the experimental uncertainty which starts 
now to be dominant when estimating the significance for Aa^. 
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Figure 2: The Pion Form Factor prediction based on r data and PDG information. The most 
important experimental data are superimposed; they do not influence the predicted curve. 
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Figure 13: The deviation Aa^ x 10 between experiment and theory for the anomalous mag- 
netic moment of the muon. The leading hadronic vacuum polarization contribution has been 
estimated via global BHLS model fits for different e+e^ — )• tt+tt^ data samples. The r predic- 
tions - using or not the spacelike data - are given in the top pair lines, followed by the fit results 
using each scan or ISR data sample in isolation or combined; see text for comments. 
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Figure 14: A set of recent estimates of the muon anomalous magnetic moment deviation from 
the BNL average value [f39l l40ll . Our own (updated) estimates are figured by BHLS:: A and 
BHLS::B for respectively configurations A and B. The results obtained using all the ISR data 
samples inside the global fit procedure are displayed under the tag Full BHLS. The statistical 
significance of each Aa^ is displayed on the right hand side of the Figure. See the text for more 
information. 



[35] NA7, S. R. Amendolia etal, Nucl. Phys. B277, 168 (1986), A Measurement of the Space 
- Like Pion Electromagnetic Form-Factor. 

[36] E. B. Dally et at., Phys. Rev. Lett. 48, 375 (1982), Elastic Scattering Measurement of the 
Negative Pion Radius. 

[37] ATLAS, G. Aad et al., Phys. Lett. B716, 1 (2012), 1207.7214, Observation of a new 
particle in the search for the Standard Model Higgs boson with the ATLAS detector at the 
LHC. 

[38] CMS, S. Chatrchyan et al., Phys. Lett. B716, 30 (2012), 1207.7235, Observation of a 
new boson at a mass of 125 GeV with the CMS experiment at the LHC. 

[39] Muon G-2, G. W. Bennett et al., Phys. Rev. D73, 072003 (2006), hep-ex/0602035. Final 
report of the muon E821 anomalous magnetic moment measurement at BNL. 



43 




- (b) 
- 




i 


\ " BaBar 












- • NA7 




7 










-0.4 



-0.4 -0.2 0.2 

s (GeV') 




-0.3 



s (GeV') 



Figure 15: The Pion Form Factor around s = 0. The tt+tt annihilation data submitted to the 
global fit are only the NSK and KLOEIO samples. The full curve is the fit function derived 
when including also the spacelike data from NA7 [i35il and Fermilab fl36l1 : the dashed curve 
is associated with the fit excluding the spacelike data. The BaBar data - shown in (b) for 
illustration - are not submitted to the fit. KLOE08 has no data point located within the plotted 
window. For clarity, only the data from NA7 [|35l are shown in (a) and (b). In (c) One magnifies 
the spacelike region and plot the data [[35l[36ll . used or not within the fit procedure. The vertical 
lines in (a) and (b) show the locations s = and s = 4m^. 
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